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Abstract 

The measurement of the inelastic photoproduction of charmonium at 
HERA seems to have ignited a new charmonium crisis. The (already discred- 
ited) color singlet model fits the data for large charmonium energy fraction 
z, where the NRQCD model qualitatively fails. We here point out that by 
the straightforward inclusion of color singlet and octet processes in the soft 
color (color evaporation) scheme, the HERA data can be accommodated for 
all z. We anticipate that the color singlet model will fail at low z, as it does 
in hadroproduction. 



I. INTRODUCTION 



The once conventional treatment of color in perturbative QCD calculations, i.e., the color 
singlet model (CSM), has run into a serious impasse in describing the data on the production 
of charmonium states at the Fermilab Tevatron J^j . Specific proposals to solve this problem 
agree on the basic elements: i) onium production is a two-step process where a heavy quark 
pair is produced first, followed by the nonperturbative formation of the asymptotic states, 
and ii) color octet as well as singlet cc states contribute to the production of charmonia, as 
clearly demonstrated by the data |],@]. Neglecting octet contributions is the source of the 
failure of the CSM. Two formalisms have been proposed to remedy this mistake. The basic 
assumption of the soft color approach HU is that no observable dynamics is associated with 
the soft processes that connect the color of the perturbative charm pair with the colorless 
charmonium bound state. This scheme, although far more restrictive than other proposals, 
successfully accommodates all features of charmonium and bottomonium production. The 
original QCD computations || referred to it as the color evaporation model (CEM), and 
identified it as the most logical way to incorporate color constraints into perturbative QCD 
calculations. The data suggests it still is. It correctly predicts the energy and final state 
momentum dependence of charmonium and bottomonium hadro- and photoproduction at 
all energies, as well as their production in electron-positron colliders. It makes the sweeping 
prediction that ^'s are produced unpolarized, a signature of the absence of dynamical final 
state effects. If such effects were discovered, there is a possible scheme, successful in the 
theoretical treatment of decays, to accommodate them: the non-relativistic QCD (NRQCD 
or COM, for color octet model) [JJ. Applying NRQCD to the bulk of the data has been a 
challenge and has, in fact, run into difficulties with recent HERA data. 

The measurement || of the inelastic photoproduction of charmonium at high energy 
has triggered a new charmonium crisis. The CSM, which qualitatively fails to describe the 
hadroproduction of charmonium, fits the HERA data for large charmonium energy fraction 
z 0], where the NRQCD model qualitatively fails ||. In this work we show that the CEM 
accommodates the data at all z. Moreover, we anticipate that the CSM will fail at low z, as 
it does in hadroproduction, since the J j ip photoproduction is mediated by resolved photons 
in that region. 

The evidence is compelling that Nature operates according to the color evaporation 
scheme. The CEM formalism simply predicts that, up to color and normalization factors, 
the energy, xf- and pr-dependence of the cross section are identical for the production of 
onium states and DD pairs. This is indeed the case @,(|,|9|]. The sum of the cross sections 
of all onium and open charm states is described by @] 



where the cross section for producing heavy quarks, a c£ , is computed perturbatively, irre- 
spective of the color of the cc pair. M c5 is the invariant mass of the cc pair. The CEM 
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assumes a factorization of the production of the cc pair, which is perturbative and process- 
dependent, and the materialization of this pair into a charmonium state by a mechanism 
that is nonperturbative and process independent. This assumption is reasonable, given that 
the characteristic time scales of the two processes are very different: the time scale for the 
production of the pair is the inverse of the heavy-quark mass, while the formation of the 
bound state is associated to the time scale 1/Aqcd- This approach to color is also suggestive 
of the unorthodox prescription for the production of rapidity gaps in deep inelastic scattering 
and in the production of hard jets at the Tevatron [p!0| -[T^|. 

Comparison with the ip data requires knowledge of the fraction of produced onium 
states that materialize as ip's, i.e., 

&ifi Pif) ionium ) (3) 

where p^ is assumed to be a constant. This assumption is in agreement with the low-energy 



data [0,|J] • Notice that a single nonperturbative factor p describes a given charmonium 
state, regardless of the spin and orbital angular momentum of the charm pair. Analysis of 
charmonium photoproduction allowed us to determined that p^ ~ 0.43-0.5 0. Even this 
parameter can be accounted for on the basis of statistical counting of final states JT5| ], as 
expected in a scheme with no final state dynamics. 

Quantitative tests of color evaporation are made possible by the fact that all ^-production 
data, i.e. photo-, hadropro duct ion, Z-decay, etc., are described in terms of the single pa- 
rameter p^ describing the frequency by which a charm pair turns into a J/ip via the final 
state color fluctuations. Its value should, in the spirit of the model, be the same for all pro- 
cesses. We have demonstrated the quantitative precision of the color evaporation scheme 
by showing how it accommodates all measurements, including the high energy Tevatron 
and HERA data, which have represented a considerable challenge for other computational 
schemes. Its parameter-free prediction of the rate for Z-boson decay into ^'s is an order of 



magnitude larger than the color singlet model and consistent with data [16 . 



II. RESULTS AND CONCLUSIONS 

Application of the the soft color scheme to inelastic charmonium photoproduction is 
straightforward. The photoproduction cross section is related to the ep cross section by 



T ep^J/*l>X\ 
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where (t 7P -^j/^x(W) is a function of the center-of-mass energy of the 7p system, and <E> is 
the distribution function of photons with virtuality Q 2 . The Q 2 integration ranges from 
Qmin = M lv 2 /{ 1 ~ V) to Qliax ( = 4 GeV 2 ). The limits for the y integral are y max / min = 
max i mini ' s i "where W max / min represent the maximum/minimum experimental values of W 
in the range 40 < W < 140 GeV 0. We used the following photon distribution function: 
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According to the parton model, the cross section for the J / ip photoproduction at a given 
center-of-mass energy W is 



v-yp->JMx(W) = J dx A J dx B fAh{x A )f B/p {x B ) &ab->j/tPx{s) 
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where the subprocess cross section a is given by Eqs. ([!]) and ([3]) in the CEM. Here, \f§ = 
■sJxaXb W is the center-of-mass energy of the subprocess AB — - > J/ipX. /^/ 7 (fs/p) is 
the distribution function of the parton A (B) in the photon (proton). For direct photon 
interactions (A = 7) we have f A ^(xA) = 5(xa — 1). 

The fraction of photon energy transferred to the J/tp in the proton rest frame is given 
by the Lorentz invariant: 

* = (7) 
7 v 

where Pj/^ i7iP is the four-momentum of the J/ip, photon, and proton, respectively. Following 
the cuts applied by the experiments 0, we impose that the transverse momentum of the 
J/ijj is larger than 1 GeV. 

jg — ► cc is the leading order process. It is, however, only important for large values of z. 
It was, in fact, the process used in Ref. to empirically determine the value of p^. For the 
range of z-values we are focusing on here, the direct photon contribution is dominated by 

7 g -> g c c , 
7 q q c c , 

see Fig. ^. The charm quark pair in jg fusion is produced in color singlet and octet config- 
urations, while 75 fusion leads to colored pairs only. These processes have to be evaluated 
with some caution because for soft values of the gluon transfer momentum t, see Fig. [1], the 
diagram represents the QCD evolution of the initial state distribution functions. In order 
to avoid double-counting one imposes \t\ < (2m c ) 2 ||17|| . We also included resolved photon 



processes, which proceed via quark-quark, quark-gluon, and gluon-gluon fusion into cc+ 
quark (gluon). The scattering amplitudes were evaluated using the packages MADGRAPH 



TJ and HELAS [I| 



The resulting charmonium photoproduction cross section is shown in Fig. y as a function 
of the center-of-mass energy W for 0.4 < z < 0.8. We fixed = 0.43, and used GRV 
distribution functions for the proton [2(J and photon [^T|. The adopted renormalization 



scale is the charm-quark mass m c for A 1 - 4 ) = 300 MeV, while the factorization scale is V^- 
We varied m c between 1.2 and 1.4 GeV. The relatively large uncertainty in our prediction 
indicates the need to include higher order corrections. This is expected because 0{a em a 2 s ) 
processes, although formally higher order, represent the dominant mechanism for z / 1 
charmonium production. Nevertheless, it is important to notice that it is not a challenge to 
describe the data. 

Having fixed any freedom in the definition of the perturbative series, we confront the 
CEM calculation with the ^-dependence of the cross section for the inelastic J / ip photopro- 
duction in Fig. H This figure shows that CEM successfully passes the test where NRQCD is 
at variance with the data M. The result supports the underlying assumption of CEM that 
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the probability of charmonium formation from perturbative cc pairs is independent of the 
color state, angular momentum and spin of the charm quark pair. 

Finally, we show in Fig. |4] the direct-photon, resolved-photon, and total contributions to 
the z-dependence of the inelastic J/ ip photoproduction. The inclusion of resolved processes 
strongly increases the J/i/j yield when z — > 0. The measurement of the J/ip production in 
this region can seal, once more, the fate of the CSM, since we anticipate that the omission 
of color octet gluon-gluon fusion processes will result in too small cross sections, as already 
observed in the Tevatron data. In principle, this is remedied in the NRQCD framework 
which, however, fails for larger values of z's. 
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FIG. 1. 0{a em a 2 s ) processes contributing to direct (a) and resolved (b) charmonium photopro- 
duction. 
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FIG. 2. Inelastic J/ip cross section for p T > 1 GeV and 0.4 < z < 0.8. We used the GRV94-LO 
structure function, A^ 4 ) = 300 MeV, renormalization scale fiR = m c , and factorization scale Q 2 = s. 
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FIG. 3. Inelastic J/tp production as function of the inelasticity parameter z for pj- > 1 GeV. 
We chose the parameters as in Fig. g. 



9 



N 

D 



10 



10 



10 



Total 
Direct 
Resolved 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 4. Inelastic J/tp production as function of the inelasticity parameter z for pj- > 1 GeV. 
We chose the parameters as in Fig. |2| and fixed m c = 1.3 GeV. 
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